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Abstract. The motion of a dark soliton is investigated in a one-dimensional 
dilute Bose-Einstein condensate confined in a harmonic trap and an optical lattice. 
The harmonic trap induces a dynamical instability of the soliton, culminating in 
sound emission. The presence of the optical lattice enhances the instability, and 
in addition, dephases the emitted sound waves, thus preventing stabilisation of 
the soliton by sound reabsorption. This instability can be probed experimentally 
by monitoring the soliton oscillations under various lattice configurations, which 
can be realised by changing the intensity and angle between the laser beams 
that form the lattice. For short enough times, such that the emitted sound does 
not reinteract with the soliton, the power emitted by the soliton is found to 
be proportional to the square of the local soliton acceleration, which is in turn 
proportional to the deformation of the soliton profile. 
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1. Introduction 

The realisation of dilute atomic Bose-Einstein condensates (BEC) in multi-dimensional 
optical lattices has opened the door to study novel quantum effects within a pure and 
highly tunable system. These include the quantum phase transition from superfluid 
state to Mott insulator [I], interference of coherent matter waves 0, Bloch oscillations 
[3], number squeezed states |3J, and arrays of Josephson junctions 0. Following the 
recent experimental creation of dark solitary waves in harmonically trapped BEC 
HO El IE] it is a now a realistic prospect to consider these macroscopic excitations in an 
optical lattice geometry. 

A dark soliton is a one-dimensional localised wavepacket consisting of a notch 
in the ambient density and a phase slip across its centre. It is supported in 
repulsive/defocussing nonlinear media where the shape of the wavepacket is preserved 
by a balance between dispersive effects and nonlinearity. Dark solitons have been 
investigated thoroughly, both experimentally and theoretically, within the analogous 
field of nonlinear optics [§]. However, the nature of trapped BEC poses new questions 
for the dynamics and stability of dark solitons. Thermal instabilities, arising from 
interaction with the thermal cloud, dissipate energy from the soliton, causing it to 
accelerate to the edge of the condensate and ultimately disappear J0|> while slow 
solitons are sensitive to quantum fluctuations 

A dark soliton experiences two kinds of dynamical instabilites. Firstly, a dark 
soliton, which is strictly a one-dimensional object, is prone to transverse instabilities 
[121 1131 114) when embedded in a two or three-dimensional system. For example, 
dark solitary waves under weak transverse confinement have been observed to bend 
via the "snake" instability and ultimately decay into vortices, in both BECs jTH] and 
nonlinear optics . However, by employing tight transverse confinement to freeze out 
the transverse degrees of freedom, quasi-one-dimensional condensates can be created 
|17lll8irHni20| . in which dark solitons are expected to have greatly enhanced lifetimes 
|21| . The second dynamical instability arises from the fact that a dark soliton is a 
solution of a homogeneous system. The longitudinal confinement featured in BEC 
experiments breaks the integrability of the system, rendering the soliton unstable 
to motion through the inhomogeneous background density, and leading to decay via 
sound emission |22l I23| . In previous work we studied this decay mechanism within a 
harmonic geometry, and isolated two limiting cases. If the emitted sound can escape 
the system, the power radiated by the dark soliton is proportional to the local soliton 
acceleration squared |24| . which in turn is intimately linked to the deformation of the 
soliton profile caused by the emitted sound waves [23]. In the opposite case of an 
infinite harmonic system, sound reabsorption occurs and the soliton decay becomes 
fully stabilised. To control this experimentally, we suggested creating a soliton in a 
tight dimple trap within a weaker harmonic trap. In such a geometry, the depth of the 
dimple (relative to the chemical potential) provides a sensitive handle on the soliton- 
sound interactions, which can be observed experimentally by monitoring changes in 
the soliton amplitude and period |24j . 

Recent work by Frantzeskakis et al. |26H27| has discussed the dynamics of slowly- 
moving dark solitons within optical lattices, including the effects of sound emission. 
The entire range of lattice periodicities were considered in [2Z|, and, where appropriate, 
analytic models for the soliton dynamics were derived. In our work we perform detailed 
quantitative measurements of the sound emission from the soliton and show that the 
acceleration squared law discussed above is valid (prior to reabsorption of sound). We 
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concentrate on much faster solitons than |27| . and discuss in detail a range of lattice 
periodicities that can be readily achieved experimentally, where the lattice wavelength 
is comparable to, or slightly larger than, the soliton width. Given the number of 
current BEC experiments featuring optical lattices |21 El EI 03 123 and the rapid 
soliton decay that can be induced in certain limits, this system may be preferential to 
investigate soliton-sound interactions than the dimple trap proposed in |24| . 

In Section 2 we introduce the system under consideration, the basic properties of a 
dark soliton, and the Gross-Pitaevskii equation used to simulate the soliton dynamics. 
In Section 3 we illustrate the dynamical behaviour of a dark soliton in an optical 
lattice, in particular the emission of sound, and compare this to the behaviour in the 
absence of the optical lattice. The dependence of the soliton dynamics on the lattice 
height and periodicity is discussed in Sections 4 and 5, and in Section 6 we analyse in 
detail the soliton dynamics and quantify the power radiated. Finally, in Section 7, we 
make some concluding remarks. 



2. Theory 

The dynamics of a trapped atomic BEC at low temperatures can be well described 
by the mean-field Gross-Pitaevskii (CP) equation. Assuming a quasi-one-dimensional 
condensate (oriented along the x-axis), the GP equation reduces to its one-dimensional 
form, 

^a^M) = (_tf_cP_ + y [x) + gmXt t)? _ ^ ^ t y (1) 

where ip( x ) is the mean-field order parameter, m is the atomic mass, and /i = ng is the 
chemical potential of the effectively one-dimensional system. The nonlinear coefficient 
g arises from the interatomic interactions and has an effective one-dimensional form, 
g = 2h 2 a/(ml 2 J, where l±_ = \fn~j (muj±) is the transverse harmonic oscillator length 
and a is the three-dimensional s-wave scattering length. Repulsive interactions, g > 0, 
are generally necessary to support a dark soliton, and will be assumed from now on. 
However, dark soliton-like structures can, under appropriate conditions, be supported 
by attractive interactions |2H| . 

We consider an external potential of the form, 

V(x) = ^mojlx 2 + V sin 2 (27re/A). (2) 

The first term represents a harmonic magnetic trap with frequency lu x and the 
second term denotes an optical lattice with height Vq and periodicity A/2, where 
A is determined by the wavelength and relative angle between the laser beams forming 
the lattice. Figure 1 shows a schematic of the trap configuration and initial condensate 
density. 

On a uniform background n , the wavefunction of a dark soliton with speed v 
and position [x — vt) is given by, 

, (x — vt) 



ip(x,t) = y/n^er 1 ^^ [ /3tanh 



ft- 



(3) 



where j3 — yl — (v/c) 2 and £ = is the condensate healing length 

corresponding roughly to the size of the soliton PHj. The soliton speed v depends 
on its depth rid relative to the background density and the phase slip S across its 
centre via v = y/n^ — rid = ccos(S/2), with the maximum speed set by the Bogoliubov 




Figure 1. Schematic of the initial condensate density (solid bold line), in relation 
to the peak condensate density no , under the confinement of a harmonic trap and 
optical lattice (dashed line), with a non-stationary soliton located at the origin. 
The optical lattice has height Vq and wavelength A. The corresponding density 
in the absence of the optical lattice is also illustrated (solid grey line). 



speed of sound c = */ /x/m. The renormalised soliton energy (after subtracting off the 
background fluid contribution) is given by UJ, 




A stationary dark soliton reaches zero density (n<2 = no), features an abrupt 7r phase 
slip, and maximum energy, whereas a soliton with maximum speed c has no phase or 
density contrast, and therefore zero energy. Note that, to first order, a dark soliton 
behaves dynamically like a classical particle |31l I32| . 

Throughout this paper we employ harmonic oscillator units, where length is 
measured in terms of the harmonic oscillator length l x = \Jfij (mwj and time in 
units of . In addition a fixed harmonic trap is always considered, corresponding 
to jj, = 70Hu) x , and simulations are performed for fixed chemical potential. We 
will assume realistic experimental parameters: lu x — 2tt x 5Hz, lo±_ = 250^, 
and 113D = 8huJ, where uj = (cj^lj^) 1 / 3 . Then, for a 87 Rb ( 23 Na) condensate, 
our space and time units correspond to l x = 4.8 (9.4) /im and uo^ 1 = 0.03ms. 
Furthermore these parameters roughly relate to a healing length £ = 0.58 (1.13) /im, 
speed of sound c = 2.5 (1.3) mm/s, one-dimensional peak condensate density 
n 1D = 2 x 10 7 (5 x 10 7 ) m _ \ and total number of atoms N = 4000 (18000). All 
effects considered in this paper manifest themselves on timescales smaller than the Is 
anticipated soliton thermodynamic lifetime |10|. 

In order to realistically probe a range of optical lattice periodicities, we envisage 
employing two almost co-propagating laser beams, of wavelength Ai asor . By varying the 
angle between the beams 9, the effective lattice wavelength can be increased beyond the 
actual laser wavelength, A — * Ai as0 r/ sinf?. Assuming a laser wavelength Ai aS er ~ 0.8/im 
(typically produced by a diode laser), then lattice wavelengths up to the order of 
A ~ 20^m can easily be produced, corresponding to A ~ 5l x for 87 Rb. 

3. Soliton dynamics in an optical lattice and sound emission 

We first consider the behaviour of a dark soliton under harmonic confinement only, 
in order to compare and contrast this to the behaviour observed when an optical 
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Figure 2. Space-time carpet plots of the renormalised condensate density (actual 
density minus background density) for (a) harmonic confinement only and (b) the 
same harmonic trap plus an optical lattice with Vb = 0.2/Li and A = 1Al x . Both 
cases start with a dark soliton at the origin with speed v = 0.5c. Sound waves 
are clearly emitted by the dynamically unstable soliton, not only as it climbs the 
trap walls but also as it traverses the lattice sites. Note that the sound waves in 
(b) have a much greater amplitude (see the different density scales). 



lattice is additionally imposed. Figure 2(a) shows a space-time carpet plot of the 
renormalised condensate density (actual density minus background density) of a dark 
soliton, initially at the origin with speed v = 0.5c. The soliton oscillates back and 
forth in the trap at a rate roughly equal to u s = u) x I y/2, as expected [El E2 E3 ESj • 
As it climbs the side of the trap, passing through the inhomogeneous background 
density, it emits counter-propagating sound waves of opposing amplitude, and becomes 
asymmetrically deformed [35] . The sound reflects off the edge of the condensate and 
subsequently reinteracts with the soliton around t ~ 5W" 1 . An equilibrium is rapidly 
established where the soliton, on average, reabsorbs as much sound as it emits. Hence 
the soliton energy features oscillations due to the emission/reabsorption process, but 
with no net decay [33|. The soliton position and energy as a function of time are 
shown respectively by grey lines in figure 3(a) and (b). 

Consider now the case of imposing, in addition to the harmonic confinement, an 
optical lattice with height Vq — 0.2/x and wavelength A = 2Al x , as shown in figure 2(b). 
The presence of the optical lattice causes the soliton dynamics to differ significantly 
from the pure harmonic case. As the soliton oscillates, it not only climbs the walls 
of the harmonic trap, but also traverses the lattice bumps. It thereby experiences 
increased dynamical instability and emits additional sound waves as it passes over 
each lattice site. The harmonic trap still contributes to the sound emission, as seen 
by the faint 'white' band of sound in figure 2(b), but this emission is heavily obscured 
(compared to figure 2(a)) by the lattice-generated sound waves. Importantly, in this 
case, reflection of sound from the edge of the condensate no longer leads to the 
establishment of an equilibrium between sound emission and reabsortion. This is 
because the presence of the optical lattice dephases, or mixes, the emitted sound such 
that only a proportion is reabsorbed. The soliton becomes faster and shallower, and 
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Figure 3. Effect of lattice height on the soliton dynamics, (a) Path and (b) 
energy of the soliton in a harmonic trap (solid grey lines), and a harmonic trap 
plus an optical lattice with A = 2Al x and lattice heights Vb = 0.05/j (long-dashed 
line), Vb = 0.1/x (dashed line), Vb = 0.2/j (dotted line) and Vb = 0.4/x (solid line). 

'anti-damps' to larger amplitudes. By t ~ the soliton has becomes so fast and 

shallow that it is virtually indistinguishable from, and gradually merges into, the sound 
field. The soliton position and energy as a function of time are shown respectively by 
dotted lines in figure 3(a) and (b). 

4. Effect of varying the lattice height 

We now consider the effect on the soliton dynamics of varying the height of the optical 
lattice Vq. Throughout this section, the periodicity of the lattice is kept constant at 
A/2 = l.2l x . Figure 3 shows the soliton position and energy as a function of time 
for various lattice heights. The soliton energy E s is calculated by integrating the GP 
energy functional, 

e (vo = £ivvf +vivf + §hAi 4 , (5) 

across the 'soliton region', conveniently defined to be (x s ± 5£), and subtracting from 
this the corresponding contribution of the background fluid (obtained from the time- 
independent solution of equation (1) in the absence of a soliton). 

Figure 3 shows soliton oscillations and energies for various lattice heights. 
Increasing the lattice height leads to two distinct effects which modify the soliton 
oscillations in opposing manners. Firstly, the optical lattice modifies the condensate 
profile, periodically reducing the density. This inhibits the motion of the soliton 
and tends to shift its turning point towards the centre of the trap and increase its 
oscillation frequency. This effect can be seen by the dashed line (Vb = 0.1//) in 
figure 3(a). Indeed, if the lattice height is sufficiently large, the soliton will initially 
be trapped within a single lattice site. Secondly, the optical lattice induces sound 
emission from the soliton. The soliton energy decays essentially monotonically, with 
the average decay rate increasing with lattice height. This decay leads to a gradual 
decrease in the oscillation frequency and increase in the amplitude. For small lattice 
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Figure 4. Effect of small lattice periodicity, A/2 ~ £ (or less), compared to no 
optical lattice, (a) Path and (b) energy of a dark soliton in a tight optical lattice, 
A = 0.3l x (2.5§), with height Vb = 0.1/x (dashed line), compared to the absence 
of the optical lattice (grey line), (c) The density (grey line) for this tight lattice 
cannot fully heal, shifting the average density down (bold line) in relation to the 
harmonic trap density (dashed line), (d) Since the fluid cannot heal fully to the 
tight lattice (grey line), the effective potential tends towards a harmonic trap 
shifted upwards (solid line) by the mean lattice potential Vq/2. The harmonic 
potential is also shown (dashed line). 



heights, e.g. Vb = 0.05/z (long-dashed line in figure 3(a) and (b)) the optical lattice 
poses little resistance to the soliton motion, and the low level of sound emission it 
induces causes the main observable effects. These are a small decrease in the soliton 
oscillation frequency and a slow decay of the soliton energy. The competing effects are 
most clearly observed for Vb = 0.4/1 (solid line in figure 3(a)). The soliton is initially 
confined to a single lattice site, with a high oscillation frequency. Sound emission 
causes the soliton to gradually 'anti-damp' out of this site into adjacent sites, leading 
to a decrease in its oscillation frequency and an increase in its amplitude. This effect 
has already been reported in |57] . The soliton rapidly accelerates towards the edge of 
the cloud where it ultimately merges into the sound field. However, in the limit Vq > /i, 
an array of mini-condensates is formed. Each lattice site becomes an effectively closed 
system, and the emitted sound cannot escape the site. In this limit the GP equation 
predicts that full sound reabsorption by the soliton occurs, which stabilises it against 
decay, as already observed in [21]. In the limit of very small lattice periodicities, the 
discrete nonlinear Schrodinger equation may be a more appropriate description, in 
which case additional oscillatory instabilities come into play |34l |2Ej ■ 

5. Effect of varying the lattice periodicity 

We next consider how the periodicity of the lattice affects the soliton dynamics, an 
effect also considered analytically and numerically in j27J. We start by discussing an 
anomalous regime, the tight lattice regime, where the lattice periodicity is of the order 
of, or less than, the soliton width. The path and energy of a dark soliton with initial 
speed v — 0.5c under the influence of a tight optical lattice (A = 0.3/a; (2.5£) and 
Vb = 0.1//) are shown by the dashed line in figure 4(a) and 4(b) respectively. The 
path of the soliton is identical to that in the absence of the optical lattice (solid grey 
line), while the energy has the same form but is shifted to a lower value. This effect 




Figure 5. Effect of lattice periodicity on the soliton dynamics, (a) Path and (b) 
energy, rescaled by the initial soliton energy E' s n,t , of a soliton in a harmonic trap 
(grey line) and an additional optical lattice with height Vq = O.lfi and various 
wavelengths: A = 0.6l x (5£) (solid line), A = 1.2?^ (10£) (long-dashed), and 
A = 4.81 x (40§) (dotted line), (c) Initial soliton energy El nlt as a function of 
lattice wavelength, compared to the Vq = harmonic case (horizontal grey line), 
(d) Loss of soliton energy AE B after one full oscillation, as a function of lattice 
wavelength. 



can be explained by careful consideration of the density perturbations caused by the 
optical lattice. In the limit where the lattice wavelength is of the order of, or less 
than, the soliton width, the fluid cannot heal to the optical lattice, and the associated 
density modulations smear together. Then the background density becomes reduced 
across the whole condensate, and in particular the peak density, no — > n' . This is due 
to the fact that, in our simulations, the chemical potential \i is fixed. If, however, the 
peak density/number of particles were kept fixed, then the chemical potential would 
increase accordingly. The soliton essentially sees a constant effective optical lattice 
potential of height Vb/2, and overall it sees the ambient harmonic potential shifted 
upwards by this amount, an effect illustrated in figure 4(c) and 4(d). In both cases of a 
pure harmonic trap and a harmonic trap modified by a tight optical lattice, the initial 
soliton speed has been set to the same value v — 0.5c (i.e. the dependence on peak 
density is removed), ensuring that the paths are identical. If we assume a Thomas- 
Fermi relationship between potential and density, n(x) = no — V(x), then we would 
expect, to first order, the peak density to be decreased by (Vo/2)no = 0.05n . Equation 
4 implies that the initial soliton energy will be rescaled by a factor of (n /no)2, i.e. 
£j, nlt — > 0.926-E""*, and this is indeed the modified soliton energy observed. 

We now investigate the effect of lattice periodicity on the soliton dynamics. The 
soliton path and energy for various lattice wavelengths are shown in figure 5(a) and 
(b) respectively. In order to remove the effect of reduced peak density for sufficiently 
tight lattices (based on our fixed chemical potential simulations), the soliton energies 
in figure 5(b) have been rescaled in terms of the initial soliton energy. This effect is 
quantified in figure 5(c) which shows the variation in the initial soliton energy as a 
function of lattice wavelength. 
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Looking now at the decay of the (normalised) soliton energy as a function of time 
for various lattice wavelengths, we find that increasing the lattice wavelength (but still 
remaining far from the tight lattice regime considered earlier) leads to the decay of the 
soliton energy and an associated increase in the soliton amplitude and frequency (e.g. 
A = 0.6Z X , solid line in figure 5(a) and(b)) in relation to the pure harmonic case (grey 
line). The rate of this decay increases up to some resonant wavelength (A = 1.22 x , 
dashed line). Any further increase in the wavelength leads to a slowing of the decay. 
This is illustrated in figure 5(b) by the reduced decay rate for A = 4.8£ (dotted line) 
compared to A = 1.2£ (dashed line). 

The dependence of the soliton decay on the lattice periodicity can be explained 
in terms of the density inhomogcneity that the soliton experiences. To illustrate this 
analysis more quantitatively, figure 5(d) plots the loss in soliton energy after one full 
oscillation as a function of lattice wavelength. For tight lattices, where the lattice 
periodicity is of the order of, or less than, the soliton width, the soliton effectively 
sees a displaced harmonic trap, and so experiences no net decay, due to complete 
sound reabsorption. In this limit, the soliton energy oscillates periodically due to the 
continuous sound emission/reabsorption process. This takes the form a beating effect 
caused by the different oscillation frequencies of the soliton and sound (dipole mode) 
|24| . Therefore, one full oscillation of the soliton in the trap does not correspond 
to one beat period, and hence the soliton energy has not yet returned to its initial 
value, explaining the small apparent energy loss. For progressively weaker lattices, the 
density begins to heal to the optical lattice, thereby modifying the condensate density 
from the harmonic trap profile. This increases the density inhomogeneity that the 
soliton sees, and induces an increase in the rate of decay. Furthermore, the presence 
of these density modulations dephases the sound field such that the soliton is no longer 
fully stabilised by sound reabsorption. 

Increasing the lattice periodicity further leads to some critical wavelength the 
condensate density can heal fully to the lattice potential. This represents the system 
with the maximum density inhomogeneity (for a fixed lattice height) and therefore 
features the most rapid soliton decay. This point is indicated by the resonance in 
figure 5(d) and occurs for A ~ 1.2l x (10£). Further increase in the wavelength does 
not change the depth of the density modulations, but merely stretches them out. This 
now reduces the density gradient experienced by the soliton [2SJ and causes a reduction 
in the decay rate. As A — > oo the decay will asymptotically tend towards zero as the 
trap potential returns to becoming effectively harmonic. 

6. Quantifying the sound emission 

We will now quantatively discuss the soliton dynamics up until the time when the 
sound reinteracts with the soliton. Figure 6 shows the short-term dynamics of a dark 
soliton in the optical lattice geometry considered in Section 3 (Vb = 0.2/U, A = 2Al Xl 
dotted line of figure 3), along with the corresponding dynamics in the absence of 
the optical lattice (grey line of figure 3). The soliton motion (figure 6(a)) is clearly 
inhibited by the modified local density (figure 6(b)), causing it to turn around earlier 
and at a lower amplitude than in the absence of the optical lattice. The soliton 
accelerates and decelerates periodically as it traverses the lattice sites, as shown in 
figure 6(c) (solid line, left axis) with the oscillation weakly modulated by the additional 
acceleration induced by the ambient harmonic trap (grey line in same figure). This 
contribution is zero when the soliton is at the centre of the trap and a maximum when 
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Figure 6. Comparison of dynamics of a soliton with initial speed 0.5c in a 
harmonic trap plus optical lattice with Vq = 0.2fi and A = 20£ (black line) and 
corresponding dynamics only in the presence of the same harmonic trap (grey 
lines), (a) Path of soliton versus time, (b) Background density at the soliton 
position, (c) Soliton acceleration (bold line, left axis) and deformation parameter 
(a; s — x cul ) (dashed line, right axis), (d) Soliton energy, (e) Power emitted by the 
soliton as computed from the GP energy functional (solid line), compared with 
the acceleration squared law of equation (7) with k = 1.8(h/l x ) (dotted line). 



the soliton is at an extremum. 

The emission of oppositely "charged" sound waves from the soliton leads to an 
asymmetric distortion of the soliton profile, and a shift of the soliton centre of mass 
x cm from the density minimum x s , where the soliton centre of mass is defined as, 

_ J>(M 2 ~n)dx 

Xcm ~ IM 2 -^ ■ 1 j 

This shift, shown in figure 6(c) (dashed line, right axis), is directly proportional to the 
soliton acceleration (solid line, left axis), in accordance with earlier findings |25j. 

The motion of the soliton through the lattice causes a periodic cascade of the 
soliton energy (figure 6(d)). This cascade is more evident when considering the power 
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emitted by the soliton, shown in figure 6(e). The power consists of a series of peaks 
of various amplitudes, which clearly occur whenever the soliton traverses the side of 
a lattice site, i.e. where the density gradient is a local maximum. The peaks are 
greatest at around t ~ \.buj~ l , which corresponds to the point of maximum soliton 
acceleration. Indeed, the power emission can be well fitted by an acceleration-squared 
law, 

where k is a constant coefficient (k = 0.02bb(h/l x ) in figure 6(e)). This is also in 
agreement with an acceleration-squared power law derived analytically for dynamically 
unstable dark solitons in the context of nonlinear optics [HS] • Although [^j deals with a 
homogeneous system in which the instabilities arise from imperfections in the medium, 
such an approach has recently been considered for dark solitons in inhomogeneous 
condensates |24II25| . When the emitted sound reinteracts with the soliton (t > Soj^ 1 ), 
it is partially reabsorbed (as evident from the negative values acquired by the power 
curve of figure 6(e)), and equation (7) ceases to be valid. For the same reasons, the 
proportionality between acceleration and the deformation parameter is no longer valid. 



7. Conclusions 



We have made a quantitative investigation of the dynamics of a dark soliton induced 
by the presence of an optical lattice. We find, in general, that the optical lattice 
significantly affects the motion and dynamical stability of the soliton. The exception 
to this is in the tight lattice regime where the fluid cannot heal to the lattice and the 
soliton effectively sees a harmonically trapped condensate. 

By considering optical lattices of varying depths and wavelengths we have 
illustrated the sensitive dependence of the soliton dynamics on these parameters, and 
explained the dependence in terms of the density inhomogeneity experienced by the 
soliton. The effect of the optical lattice, and the associated density modulations, is 
to create extra dynamical instability in the soliton, in addition to that caused by 
the ambient harmonic trap. This induces soliton decay via the emission of sound 
waves. For deep and moderately tight lattices this effect is considerable: the soliton 
rapidly decays, accelerating to the edge of the condensate, where it disappears into 
the sound field. Unlike in a pure harmonic trap, the presence of the optical lattice 
tends to modify the emitted sound such that the full stabilisation of the soliton by 
sound reabsorption is prevented. 

The soliton-sound interaction in an optical lattice geometry can be experimentally 
probed in a quasi-one-dimensional condensate, providing an alternative method to 
that suggested in [2"3| . The dynamical dissipation is observable as a change in the 
soliton amplitude and oscillation frequency, by repeated expansion of the cloud and 
subsequent absorption imaging, after a variable evolution time in the trap. The 
timescales considered in this work are consistent with the expected thermodynamic 
soliton lifetime of Is ^UJ- The intensity of the laser beams that form the optical lattice 
controls the lattice height, while varying the relative angle between the beams allows 
the full range of lattice periodicities considered here to be realised. Modification of 
these two parameters enables full experimental control of the dynamical instability of 
the soliton induced by the optical lattice. 
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